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This article presents a unified synchronization framework with application to precision 
formation flying spacecraft. Central to the proposed innovation, in applying synchroniza- 
tion to both translational and rotational dynamics in the Lagrangian form, is the use of the 
distributed stability and performance analysis tool, called contraction analysis that yields 
exact nonlinear stability proofs. The proposed decentralized tracking control law syn- 
chronizes the attitude of an arbitrary number of spacecraft into a common time-varying 
trajectory with global exponential convergence. Moreover, a decentralized translational 
tracking control law based on oscillator phase synchronization is presented, thus enabling 
coupled translational and rotational maneuvers. While the translational dynamics can be 
adequately controlled by linear control laws, the proposed method permits highly nonlinear 
systems with nonlinearly coupled inertia matrices such as the attitude dynamics of space- 
craft whose large and rapid slew maneuvers justify the nonlinear control approach. The 
proposed method integrates both the trajectory tracking and synchronization problems in 
a single control framework. 

I. Introduction 

Motivated by distributed computation and cooperation, abundant in both biological systems (e.g., fish 
swarms) and artificial machines (e.g., parallel computers), formation flying spacecraft has been a key research 
topic among many recent advancements.^ Multiple apertures flying in precise formation are expected to 
provide unprecedented image resolution, both for astronomy and reconnaissance,^ as well as unparalleled 
reconfigurability. However, many significant technical challenges must be overcome before formation fiy- 
ing interferometers can be realized (see two representative missions in Figure 1). For instance, formation 
fiight requires extensive technology development for precise attitude and position maintenance of multiple 
spacecraft. 

The objective of this paper is to introduce a unified synchronization framework that can be directly 
applied to the position and attitude synchronization, and cooperative control of formation fiight networks, 
comprised of either identical or heterogenous spacecraft. Synchronization is defined as a complete match 
of all configuration variables of each dynamical system. We also introduce phase synchronization where 
spacecraft follow an oscillatory trajectory with some phase difference between spacecraft. In particular, 
we show that we can synchronize the attitudes and positions of multiple spacecraft faster than they track 
the common position and attitude trajectories. The combined synchronization and tracking control law 
can achieve more efficient and robust performance through local interactions, especially in the presence 
of non-identical disturbances and uncertainties. Such local interactions are key to stellar formation fiight 
interferometers that depend on precision control of relative spacecraft motions, indispensable for coherent 
interferometric beam combination (see Figure 1). 

A recent review paper^'^ highlighted the three main areas for future research that have not been thor- 
oughly addressed in the spacecraft formation fiight literature: (1) rigorous stability conditions for cyclic 
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(a) NASA formation flying interferometer TPF, courtesy: (b) NASA Stellar Imager (SI), courtesy: NASA, 
NASA, http: / /planetquest.jpl. nasa.gov http: / /hires. gsfc.nasa.gov/si 



Figure 1. Representative future NASA missions using formation flying spacecraft. 



and behavioral architectures, (2) reduced algorithmic information requirements, and (3) increased robust- 
ness/autonomy. While this paper does not provide complete solutions to these three challenges, we compare 
this paper with prior work in the aforementioned areas as follows: 

1. Rigorous Stability Condition for Highly Nonlinear Time Varying Systems 

Prior work on consensus and flocking problems using graphs, particularly popular in the robotics research 
community, tends to assume very simple dynamics such as linear systems and single or double integrator 
dynamic models. Such work can be generally applied to the synchronous spacecraft position control 
problem, which fares well with linear control. In contrast, the proposed strategy in this paper primarily deals 
with complex dynamical networks consisting of highly nonlinear time- varying dynamics that are controlled 
to track a time-varying reference trajectory or leader. Examples of such a nonlinear system include the 
attitude dynamics of spacecraft for large and rapid slew maneuvers. 

Some prior works on attitude synchronization rely on tracking a common reference (leader) spacecraft 
or nearest neighbor without local information flows. ^^^^ Local coupling control laws are suggested in [15], 
but individual spacecraft are synchronized to a constant state, thus not permitting an arbitrary reference 
trajectory. This is particularly true of consensus problems on graphs.^ In [16], a new coordination architec- 
ture for both attitude and translational dynamics was proposed to incorporate leader-following, behavioral, 
and virtual-structure approaches. A decentralized control law using the virtual structure approach was pre- 
sented for a potentially large number of spacecraft in [17]. Decentralized attitude control laws with global 
asymptotic convergence were introduced in [18]. Compared with one recent work^^'^^ that studied formation 
keeping and attitude alignment for multiple spacecraft with local couplings, the present paper introduces 
both the new attitude synchronization strategy with global exponential convergence and the novel position 
consensus strategy using phase synchronization. Previously, phase synchronization of a simple planar phase 
model of fish was studied in [21]. 

It should be noted that determining stability of nonlinear time- varying dynamic network systems is more 
involved and difficult. ^^^^ Many mechanical systems exhibit nonlinear dynamics that cannot be captured 
by linearization. One might argue that most space systems are not required to follow demanding time- 
varying trajectories, thus validating linearization or linear coupling control laws. In particular, for the 
translational dynamics, a linear coupling control law can effectively stabilize the formation flying spacecraft. 
However, global exponential convergence of the attitude dynamics is achieved only through nonlinear control. 
For instance, there is increased interest in highly agile imaging spacecrafts^ that undergo wide and rapid 
slew angle changes. Moreover, in the context of nonlinear control theory, the asymptotic convergence of 
linear control, employed to stabilize nonlinear systems, may not be sufficient for demanding future mission 
requirements. In essence, ensuring exponential tracking stability for general nonlinear systems is made 
possible only through nonlinear control, and the benefit of exponential stability, in terms of improved tracking 
performance and robustness, is illustrated in this paper. 
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We introduce contraction analysis^^^^ as our main nonlinear stability tool for reducing the complexity 
and dimensionality associated with multi- agent systems, thereby deriving exact and global results with 
exponential convergence with respect to arbitrary time-varying inputs (see the Appendix for the further 
treatment of contraction theory). 

2. Reduced Information Network 

Another benefit of synchronization is its implication for model reduction. The exponential synchronization of 
multiple nonlinear dynamics allows us to reduce the dimensionality of the stability analysis of a large network. 
The model reduction aspect of synchronization, also introduced for spatially inter- connected systems in 
Ref. 2, is further generalized and strengthened in this article. This implies that once the network is proven 
to synchronize, we can regard a network as a single set of synchronized dynamics, which simplifies any 
additional stability analysis. As shall be seen later in the subsequent sections, this model reduction has to 
do with the fact that there are two time-scales associated with the coupled nonlinear dynamics. 

In addition, the proposed control laws are of a decentralized form requiring only local velocity /position 
coupling feedback for global exponential convergence, thereby facilitating implementation in real systems. 
In contrast with some previous work using all-to-all coupling or depending only on tracking the same leader 
(reference) spacecraft without local interactions, our proposed approach will not only reduce communica- 
tion burdens, but also increase the overall performance of relative formation flight through local interactions. 
Further, we mathematically prove that the synchronized maneuvers can also be achieved by sharing partial 
state information. This partial degrees-of-freedom coupling will further reduce the amount of formation state 
information needed for formation flight. Another recent work using the passive decomposition^^ described a 
strategy of decoupling into the same system representing the internal group formation shape, and the locked 
system describing the total group maneuver. 

3. Robustness Issues 

The proposed method in this paper integrates both the trajectory tracking and synchronization problems 
toward a single control framework. Although an uncoupled trajectory control law, in the absence of external 
disturbances, would achieve synchronization to a common trajectory, the proposed control strategy can 
achieve more efficient and robust performance through synchronization using local interactions, especially 
if the external disturbances vary with each spacecraft. In addition, most previous works do not discuss a 
property of robustness in detail. In contrast, we show that the proposed decentralized control law possesses 
a property of robustness to inter-spacecraft time-delays and bounded disturbances. An adaptive version of 
the proposed control law is also presented to deal with parametric uncertainties of dynamic models. 

4. Control of Relative Translational Motions 

Whereas the present paper presents the unified decentralized control method for the Lagrangian dynamics 
of both rotational and translational motions of the spacecraft, with focus on minimizing the tracking and 
synchronization errors, the literature focused only on the relative translational motions is abundant. A 
decentralized formation control law is presented in [ 31]. A disturbance accommodating control design process 
is presented in [ 32] for minimizing the total fuel consumption for the formation as well as maintaining the 
equal level of fuel consumption for each spacecraft. Ref. 33 presents a Hamiltonian approach to modelling 
relative spacecraft motion based on a derivation of canonical coordinates for the relative state-space dynamics. 
A formation control law based on Gauss's variational equations was first introduced in [34], while a recent 
paper^^ presented a new linear time- varying form of the equations of relative motions developed from Gauss's 
variational equations. A sliding mode controller for Hill's relative motion equations was presented in [36]. The 
current paper presents a synchronization control law for the nonlinear relative spacecraft dynamics, although 
the proposed approach can be applied to any Lagrangian formulation of spacecraft motions including Hill's 
equations. 

Compared with our recent paper^^ that first reported the combined synchronization and tracking control 
law for multiple robots, this paper presents the Lagrangian formulation of spacecraft dynamics (Section II), 
introduces a new strategy for oscillator phase synchronization that enables coupled rotational and transla- 
tional maneuvers (Section IV), and discusses the properties of robustness to transmission delays and distur- 
bances more thoroughly (Section V). A few examples of formation fiying spacecraft are given in Section VI 
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for validating the effectiveness of the proposed synchronization framework. 



II. Lagrangian Formulation of Formation Flying Spacecraft 

The proposed synchronization framework is devoted to the use of the Lagrangian formulation for its 
simplicity in dealing with complex systems involving multiple dynamics. We show herein that the rotational 
maneuvers of a rigid spacecraft can be written in this Lagrangian form, thereby permitting direct application 
of the proposed synchronization strategy^' to the rotational dynamics of multiple spacecraft. Without loss 
of generality, the proposed control law can be applied to the position synchronization of formation flying 
spacecraft, or more generally to the coupled translational and attitude dynamics. 

A. Lagrangian Formulation 

The equations of motion for a spacecraft with multiple degrees-of- freedom (q^ G M^) can be derived by 
exploiting the Euler-Lagrange equations: 

UiSii^^i) = Mi(qOqi - Vi{qi), — ^ ^ = n (1) 

where i, (1 < i < denotes the index of spacecraft comprising a spacecraft formation flight network, and p 
is the total number of the individual elements. Equation (1) can be written as 

Mi(qi)qi + Q(qi, qi)qi + gi(qO = n (2) 

where gi(qi) = ^^q^'^ ; and, is a generalized force or torque acting on the i-th spacecraft. 
It should be emphasized that, among many possible choices, the C matrix is defined as 

1 ^ dMij , , 1 ^ fdMik dMj 



Then, it is straightforward to show that (M^ — 2Ci) is skew-symmetric, resulting in 

x^(Mi-2Q)x = 0, l<i<p (4) 

for an arbitrary x G M^. This skew-symmetric property can be viewed as a matrix expression of energy 
conservation, which can also be explained in the context of the passivity formalism. In the remainder 
of this paper, the property in (4) is extensively exploited for stability analysis and control synthesis using 
contraction theory.^ 

We assume that the spacecraft system in (2) is fully actuated. In other words, the number of control inputs 
is equal to the dimension of their configuration manifold (n). However, we do not require the communication 
of the full state information for the purpose of synchronization, as shall be seen in Section VI-A. 

B. Attitude Dynamics of Rigid Spacecraft 

The aim of the present section is to show that we can establish a Lagrangian formulation given in (2) from the 
rotational attitude dynamics of a rigid spacecraft (n = 3). We improve the approach in Ref. 37 in two aspects. 
First, we introduce the generalized form to explore the increasing interest in agile imaging spacecraft using 
Control Moment Gyroscope (CMG)^^ and Variable Speed Control Moment Gyroscope (VSCMG).^^ Second, 
we incorporate the Modified Rodrigues Parameters (MRPs),^^'^^ in lieu of the Rodrigues parameters, to 
overcome the singularity problem at the rotation of ±180 deg. 

Using the Euler rotational equations of motion, the following equation describes the angular velocity 
vector cc; G of the spacecraft in its body axes: 

Js/c^ - (Js/c^) X (J = U -h dext (5) 

where the internal control torque u, generated either by VSCMGs or CMGs, is defined as^^ 

u = — h -h h X cj (6) 
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Based on (6), a suitable algorithm, such as the pseudo-inverse steering logic, can determine the gimbal 
angle vector 7 once the control input u is computed from the attitude control law. 

Note that the matrix Jg/c is the total moment of inertia of the spacecraft, expressed in its body frame, 
and is symmetric positive definite. Also, h and dext, all expressed in the spacecraft body- fixed frame, denote 
the total control momentum vector by CMGs, and the external disturbance torque such as the aerodynamic 
drag torque and the gravity gradient torque. We also assume that the change of Jg/c due to the CMG gimbal 
angular rate is small (i.e., Jg/c = 0). 

In the case of VSCMGs, the control momentum vector (h) and its rate (h) can be written as^^ 



h = A^I^7 -h A^I^ll -h Atl^diag(r2)7 



(7) 



where 7 and H denote the gimbal angles and wheel speeds of the VSCMGs, respectively, while Ig and are 
the moment of inertia of the gimbal structure and the wheel. Also, A^, A^, and At denote the transformation 
matrices associated with the body- frame representation.^^ 

The purpose of introducing (7) is to show that (5) may represent either the rotational dynamics of 
spacecraft with fixed-speed CMGs (constant fi) or reaction wheels (7 = 0). 

To avoid the singularity problem of the Euler angular representation, it is often preferred to use quater- 
nions to represent an angular orientation between two different coordinate frames: 



6 6 6 6 

/5i = eisin-, /52 = e2sin-, [5^ = e^sui-, (3^ = cos - 



(8) 



where e = (ei, 62, 63)^ is the Euler axis of rotation expressed in the body frame and 6 is the rotation angle 
about e. 

The modified Rodrigues parameters can be written as:^^'^^ 



etan 



q= (^1,^2,^3)^ 
Then, the attitude of the spacecraft has the following relation: 

q = Z{q)uj 

where 

,1 



6 



Z(q) 



I( 



T 

q q 



2 

(i + gi^ 



) + qq^ +S(q) 

2 



(9) 



(10) 



(11) 



q2 - qs ) 2(gig2 - qs) 

2(^2^1 + ^3) (1 - gi^ + - ^3^) 

2(^3^1 -^2) 2(^3^2 + ^1) 



and the skew-symmetric matrix function, S(x) for an arbitrary x G 



2(gig3 + q2) 
2(^2^3 - ^1) 
(1 - qi^ - ^2' + ^3'). 

^ is defined as 








-X3 


X2 


S(x) = 







-Xi 




-X2 








(12) 



Also, note that the corresponding quaternions in (8) can be obtained from the modified Rodrigues parameters 
by the following transformation 



/3i = 2gi/(l + q^q), i = l,2,3, /34 = (1 - q^q)/(l + q^q), 



while its inverse transformation can be written 

5i=A/(l+/?4), 



1,2,3 



(13) 



(14) 
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By combining (5) and (10), the following equations of motion are obtained with respect to the modified 
Rodrigues parameters for the i-th spacecraft {I < i < p) 

Mi{qi)qi + Q(qi, qi)qi = n -\- Text,i (15) 

where 

Ti = Z"^(qi)u, rext,i = Z"^(qi)dext,i 
M,(q,) = Z-^(q,)Js/c,^Z-l(q,) (16) 
= — Z~^Js/c,iZ~-^ZZ~-^ — Z~^S(Js/c,i^i)Z~-^. 

Notice that the index i has been added to Js/c,i (hence, also to and C^) to permit complex formation 
networks comprised of multiple heterogeneous spacecraft. Also, note that S(Js/c,z^z) = ^{Js/c,i'^~^{^i)^i) 
from (10). Notice that all the terms in (15-16) are left-multiplied by Z~^(q^). However, we should not cancel 
out the common term Z~^(q^) of (15-16), which would result in having a non-symmetric M^(q^). In essence, 
we established a Lagrangian formulation for the attitude dynamics of rigid spacecraft. This allows us to 
apply a wealth of nonlinear control laws to spacecraft dynamics, including the proposed control strategy in 
Ref. 22, that were originally developed for robot dynamics. As discussed in (4), the most important feature 
of (15) is to have a skew-symmetric M — 2C due to energy conservation. Indeed, we can verify that 

Mi - 2Ci = ^^Js/c,iZ-i - Z-^J,/e,,^ + 2Z-^S(Js/e,i^0Z-' (17) 

is skew-symmetric, which follows from that S(Js/c,^^i) is skew-symmetric. Without loss of generality, the 
control torque u generated by momentum wheels (7 = 0) can be defined as = — h^. Then, the S{Js/c,i^i) 
in (16) is replaced by S(Js/c,i^i + h^) in order to account for the gyro stiffening effect of the wheels. 

In the subsequent sections, the rotational dynamics formulation in (15) is used to develop a nonlinear 
synchronization tracking control law for multiple formation flying spacecraft. 

C. Relative Translational Dynamics 

If we assume that the influence of the attitude dynamics on the translational dynamics is weak and ignored, 
the translational dynamics, modeled as double integrators,^ can be easily augmented with the attitude 
dynamics in (15). Alternatively, similar to Ref. 14,42, the coupled translational and rotational motions 
of formation spacecraft can be written in the Lagrangian form in (2). Then, the proposed decentralized 
tracking control can be effectively applied without loss of generality. For arbitrary translational dynamics, 
synchronization corresponds to xi = X2 = • • • = Xp where x^, 1 < i < p connotes a vector of biased variables 
constructed from the position vector such that Ti{t) = x^(t) + hi{t) and the separation vector hi{t) is 
independent of the dynamics. Conceptually, such a method follows that each position vector x^ can be 
defined from virtually shifted origins. In this paper, the goal is to take a different approach in which the 
phase differences of the position variables can synchronize on a spiral or circular trajectory. 

In pursuit of this goal, this section presents relative translational dynamics applicable to formation flight 
on a circular or spiral configuration in low earth orbit (LEO) such that the formation translational control 
is based on the relative dynamics with respect to the desired formation center of mass (cm.). For deriving 
a relative position control law, it is more advantageous to work in the non-inertial orbital coordinates F^^ 
as illustrated in Fig. 2. The orbital frame F^^ is defined in such a way that its origin is attached to the 
center of mass of the formation with its y-axis aligned with the position vector Rq representing the position 
of the formation center of mass in the Earth Centered Inertial (ECI) frame. The z-axis points towards the 
orbital plane normal and the x-axis completes the right hand system (see Fig. 2). 

Although this derivation can also be done in an orbital frame attached to the leader satellite, the leader 
satellite will be maneuvering in general so its orbital velocity would not be a constant. Therefore, it is more 
accurate to use an orbital reference frame attached to the center of mass of the formation as opposed to the 
leader satellite. 

For simplicity, we will assume a circular reference orbit for the formation center of mass for which the 
angular velocity of the reference orbit is simply: 

(18) 
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Formation cm. 



Figure 2. Geometry of different reference frames 



where /ig is the gravitational constant of Earth (398,600.4418 x 10^ m^/s^) and Rq is the radius of the 
formation's center of mass orbit. 

Even though we consider the J2 effects as a source of disturbance for simulation studies in this paper, 
we ignore the J2 terms for the control law development. Then, the relative dynamics of i-th satellite with 
respect to satellite-k in the orbital frame F^^ can be written as^ 



xd 



Xi - 2ujQyi -ujQ Xi^ —To- = 

itf m 



(19) 



Rf m 



where F^^ Fy^ and Fz denote the control forces while F^d^ Fyd, and Fzd represent the external disturbance 
forces, all expressed in the orbital frame F^^ . In addition, the distance between the Earth center and the 
i-th satellite is defined as 

Ri - 



x1 + (yi + i?o)2 + z} (20) 

Ignoring the disturbances, (19) can be written in a Lagrangian form with a constant mass matrix, similar 
to (2) 

(21) 
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(22) 



We make two important remarks regarding (21). First, as opposed to the attitude dynamics of spacecraft 
in the preceding section, M and C are constant, which simplifies the stability proofs. Nonetheless, M — 2C 
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is skew- symmetric, which unifies our control law design for both attitude and translational dynamics. Note 
that the C matrix in (21) is not obtained from (3). Second, if the difference between Ri and Rq is reasonably 
small, D(ri) becomes a constant matrix such that D = diag(0, 0, cJq). This may further simplify the proposed 
control law in Sec. IV. 

III. Decentralized Nonlinear Control for Attitude Synchronization 



We consider the attitude synchronization of multiple spacecraft following a common time- varying trajec- 
tory in this section. 




(a) network of identical spacecraft via communication coupling (left) or relative (b) non-identical spacecraft 

sensing (right) 




(c) concurrent synchronization of multiple spacecraft groups 



Figure 3. Formation flying networks of identical or nonidentical spacecraft using local couplings considered in 
this paper: (a) and (b). They are on balanced bi-directional regular graphs, but a more complex geometry or 
leader- follower network can also be constructed by concurrent synchronization:^^ (c). Also, without loss of 
generality, the bi-directional couplings can be extended to uni-directional couplings. 
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A. Attitude Synchronization of Spacecraft 



The following decentralized tracking control law with two- way-ring symmetry is proposed for the i-th space- 
craft in the network comprised of p spacecraft (see Figure 3(a-b)): 



Ti = Mi{qi)qr,i + Q(qi, qi)qr,i + gi(qi) - Kis^ + Kss^-i + K2Si+i 



(23) 



where a positive-definite matrix Ki G M^><^ is a feedback gain for the i-th satellite, and another positive- 
definite matrix K2 G M"^><^ is a coupling gain with the adjacent members — 1 and For two-spacecraft 
networks, the last coupling term with the i + 1-th member in (23) is not used. It should be emphasized that 
the assumption of the bi-directional coupling in (23) can be relaxed without loss of generality to account for 



a regular digraph. 



Also, cir,i and are defined such that 
= qrf(t) + A(qrf(t) - q^), = q^ - q. 



qi - qd(^) + A(q^ - qd(t)) 



(24) 



where A is a positive diagonal matrix. The time- varying desired trajectory qd(t) can be the formation fiying 
guidance command or the trajectory of a leader spacecraft. 

If we assume that a relative attitude metrology system, similar to Ref. 44, is available in addition to each 
spacecraft's own attitude measurement with respect to the inertial frame, the relative attitude errors (e.g., 
q^+i — q^) can be computed. In turn, nonlinear observers can estimate the velocity errors (e.g., q^+i — q^), 
we can rewrite (23) as 



Ti =Mi{qi)qr,i + Q(qi, qi)qr,i + gi(qO - (Ki - 2X2)3^ 

-h K2 [(q^-i - Cii) + A(qi_i - q^)] -h K2 [(q^+i - q^) + A(qi+i - q^)] 



(25) 



This can be straightforwardly applied to the translational dynamics with relative distance measurements,^^ 
as shall be seen in Section IV. This implies that such a relative sensing system can eliminate the need for 
exchanging the state information of each spacecraft; no inter-satellite communication would be required. In 
essence, we can implement the proposed decentralized control law by either the communication links or the 
relative metrology system, as illustrated in Fig. 3. 

It should be noted again that one of the main contributions of this paper lies with the use of a new 
differential stability framework, yielding the exact proof of nonlinear stability under a variety of conditions, 
while many others can still come up with nonlinear control laws similar to (23), without a rigorous stability 
proof. Note that the above control law requires only the coupling feedback of the most adjacent spacecraft 
(i — 1 and for exponential convergence (see Figure 3). Consequently, the last {p-th) satellite is connected 
with the first satellite to form a ring network as suggested in Ref. 28. In order to construct a more complex 
geometry rather than a single ring network (e.g., see Fig. 3(c)), concurrent synchronization^^ can be used as 
expanded upon in Ref. 4,22. 

Without loss of generality, we can extend the proposed control law to adaptive control^^ 



Ti = Yiki - KiSi + K2Si_i + K2Si+i 

= Miq_r,i + CiCir,i + gi(qi) - KiSi + K2S^_i + K2S^+i 
The parameter estimate a^ for the i-th member is updated by the correlation integral: 

ki = -TYTsi 



(26) 



(27) 



where F is a symmetric positive definite matrix. The stability proof of the adaptive control law, presented 
in Refs. 22, does not alter the main proofs in this paper except that the convergence result reduces to 
asymptotic in stead of exponential. Hence, we will only focus on the general control law in (2) for the sake 
of simplicity. 

The closed- loop dynamics for the whole formation, by using (2) and (23), can be written as 



[M]x+[C]x+[L^^__kJx = 



where 





"Mi(qi) ••• " 




Ci(qi,qi) ••• 






[M] = 


_ ••• Mp(qp)_ 


, [C] = 


••• Cp(qp,qp)_ 


, x = 





(28) 



(29) 
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Also, the pnxpn block matrix [L^^ -K2] -^i diagonal matrix elements, neighbored by — K2. In 

other words, from the definition of the controller in (23), [Lki -K2] ^^^^ three nonzero matrix elements 
in each row (i.e., Ki, K2, K2). 

The network graphs illustrated in Figure 3 are balanced due to bi-directional coupling.^ However, it should 
be noted that the matrix [L^^ is different from the standard weighted Laplacian found in Ref. 8. By 
definition, every row sum of the Laplacian matrix is zero. Hence, the Laplacian matrix always has a zero 
eigenvalue corresponding to a right eigenvector, 1 = (1, 1, • • • , 1)^.^ In contrast, a strictly positive definite 
[L^^ -K2] required for exponential convergence for the proposed control law in this paper. 

We present the main theorems of the paper. First, the following condition should be true for exponential 
convergence to a common desired trajectory qd(t). 



Theorem III.l Global Exponential Convergence to the Desired Trajectory 

If [L^^ -K2] uniformly positive definite, then every member of the network follows the desired trajectory 
cid{t) exponentially fast regardless of initial conditions. In other words, if Ki — 2K2 > 0, then q^, (i = 
1,2, ••• ,p, p > 3) converges to qd(t) exponentially fast from any initial conditions. For two-spacecraft 
systems (p = 2), Ki — K2 > needs to be true instead. 

Proof We present a new proof that emphasizes the hierarchical combination structure of the proposed 
control law. The equation (28) of the closed-loop dynamics corresponds to a conventional tracking problem. 
Recalling Theorem VII. 2, we construct the following hierarchical virtual system of yi and y2'- 



[M] 
I 




[C] + [L^„_K2] 
-I [A] 




(30) 



where the virtual system of yi is obtained by replacing x with yi in (28), and the system of y2 is from the 
definition of the composite variable in (24). Also, [A] = diag(A, • • • , A). 
It is straightforward to verify that (30) has two particular solutions: 




and , where |q| 

\y2 = or 



The diff"erential virtual length analysis with respect to the uniformly positive definite metric. 



[M] 
al 



3a >0 



yields 
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[M] 






al 




[M] 




T r 



-2[LL,-K2] 
2al 



where we used the skew- symmetric property of [M] — 2[C]. 
The symmetric part of the matrix B is 




:(B + B^) 



"2[Lki -K2 

al 



al 

-2a[A] 



(31) 



(32) 



(33) 



(34) 
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According to Theorem VII.l, (30) is contracting if the symmetric matrix is uniformly negative definite. 
We can always find a > such that is uniformly negative definite: 



"I<4[L^,^,-kJ[A], [L^,^,_kJ>0, and A>0 



(35) 



Accordingly, all solutions of (30) converge to each other exponentially fast, resulting in global exponential 
convergence of q to q_d{t)- The positive-definiteness of [L^^ -K2] corresponds to Ki — K2 > for two- 
spacecraft systems {p = 2). For a network consisting of more than two spacecraft (p > 3), it can be shown 
that Ki — 2K2 is a sufficient condition of the positive-definiteness of [L^^ given Ki > 0, K2 > 0. 



We now focus on the synchronization of multiple spacecraft dynamics. First, we introduce an orthogonal 
matrix, whose column vectors constitute a superset of the flow- invariant subspace of synchronization^^' 



such that V^^^x = 0. 

Since [L^^ -K2] ^ constant real symmetric matrix, we can perform the spectral decomposition 



[LL,-kJ=V[D]V^ 



(36) 



where [D] is a block diagonal matrix, and the square matrix V is composed of the orthonormal eigenvectors. 
Since the symmetry of [L^^ -K2] Si^cs rise to real eigenvalues and orthonormal eigenvectors,^^ we can verify 
that V^V = VV^ = lpn, 

The modified Laplacian [L^^ -K2] defines a regular graph, where each member has the same number of 
neighbors (= 2 for p> 3). A regular graph has the block column identify matrix [1] = [ly^,!^^, • • • 
as its eigenvectors associated with the tracking convergence rate A(Ki — 2K2) for p > 3. Hence, we can 
define a x (p — l)n matrix Ysync constructed from the orthonormal eigenvectors other than [1] such that 



V^V 




sync 



I 

* syncl 



1] 



L-*-J V sync 

V 

^ sync ^ sync_ 



(p— l)nxn 



^x (p—l)n 
[(p-l)n 



(37) 



where we used the orthogonality between [1] and Vg 



The synchronization of multiple spacecraft qi = q2 = • • • = qp corresponds to 



ync 



(38) 



Theorem III. 2 Synchronization of Multiple Identical or Heterogeneous Spacecrafts^ 

A network of p spacecraft synchronize exponentially from any initial conditions if 3 diagonal matrices Ki > 0, 
K2 > such that 

^fync [Lki,-K2]^s?/^c > 

In addition, A is a positive diagonal matrix defining a stable composite variable = q^ + Aq^ with q^ = 
c[i-c[d{t)- 

If we have uni-directional couplings on a regular graph, the preceding conditions are replaced by 
^sync (^Ki -K2] ^ -K2]"^) ^sync > (scc Rcf. 22 for details). This theorem corresponds to synchro- 
nization with stable tracking. Multiple dynamics need not be identical to achieve stable synchronization. 
It should be noted that Theorem III. 2 can be used, regardless of Theorem III.l. For example, [L^^ -K2] 
might be semi-positive definite, thus yielding indifferent tracking dynamics. In this case, we do not need a 
common reference trajectory for the synchronization of multiple spacecraft, and the spacecraft synchronize 
to the average of the initial conditions (see Ref. 22 for details). It is useful to note that the above condition 
corresponds to Ki + K2 > for two-spacecraft and three-spacecraft networks (p = 2, 3). 

Note that we can make the system synchronize first, then follow the common trajectory by tuning the 
gains properly. This indicates that there exist two different time-scales in the closed-loop systems constructed 
with the proposed controllers. For two-spacecraft systems, the convergence of exponential tracking is pro- 
portional to Ki — K2 whereas the synchronization has a convergence rate of Ki -hK2. This multi-time-scale 
behavior will be exploited in the subsequent sections. 
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B. Proof of Exponential Synchronization 



We summarize the proof of Theorem III. 2 for the exponential synchronization of multiple nonlinear dynamics, 
first reported in Ref. 22. The key result can be generalized for an arbitrary number of spacecraft, even to 
more complex structures beyond a standard ring geometry. 

Suppose that M(q) is a constant inertia matrix M, thereby resulting in C(q, q) = 0. Then, we can 
easily prove Si and S2 tend to each other using Theorem VII. 4. On the other hand, the difficulties associated 
with nonlinear time- varying inertia matrices can be easily demonstrated. In essence, M(qi) ^ M(q2) 
makes this problem intractable in general. We now present a solution to this open problem, focused on the 
synchronization of multiple spacecraft with non-constant nonlinear metrics. 

Recall the closed- loop dynamics given in (28). Pre-multiplying (28) by and setting x = VV^x result 

in 



(V^[M]V) z + (V^[C]V) z + [D]z = 



(39) 



where z = V-^x. 

Then, we can develop the squared-length analysis similar to (33) 
symmetric positive definite since [M] is symmetric positive definite. 

Using (37), the block diagonal matrix [D], which represents the eigenvalues of [L^ 
tioned from (36) 



Notice that (V^[M]V) is always 
j^^], can be parti- 



^stync [I^Ki ,-K2 ] ^ sync_ 



Di 

(p— l)nxn 







nx (p— l)n 

Do 



(40) 



It should be noted that Di = Ki — 2K2 for p > 3 (or Di = Ki — K2 for p = 2) represents the tracking 
gain while D2 corresponds to the synchronization gain. We can choose the control gain matrices Ki and K2 
such that 

D2 = Vr,„,[L^^^_KjV,,„e > Di = [1]^[L^,,_kJ[1] (41) 

This will ensure that multiple spacecraft synchronize faster than they follow the common desired trajectory. 
In other words, multiple spacecraft synchronize first, then they converge to the desired trajectory while 
staying together. 

Now let us write the virtual system of y by replacing z in (39) with y. 



(V^[M]V) y + (V^[C]V) y + [D]y = 

The above system has the following particular solutions: 



(42) 




[l]^x ^ 

^ sync^J 



and y 




(43) 



Now we need to prove that the system in (42) is contracting (i.e., (5y ^ globally and exponentially) 
in order to show that those two solutions tend exponentially to each other. Performing the squared-length 
analysis with respect to the symmetric positive definite block matrix V-^[M]V as the contraction metric 
yields 




[l]^[M][l] 
[M][l] 



* sync 



[M]V,j,„c 



^ sync 





Di 






D2 




(44) 



where we used the skew-symmetric property of (V^[M]V) -2(V^[C]V). 

If [D] > 0, or equivalently [Lki,-K2] ^ ^' combined virtual system in (42) is contracting. This in 
turn implies that all solutions of y tend to a single trajectory. As a result, [l]^x = (si + • • • + Sp)/y^ and 
V^^^x tend exponentially to zero. Note that Si, • • • ,Sp 



What is new here is the proof of synchronization x 



has already been proven for Theorem III.l. 
0. Note that the synchronization occurs even 



if the tracking dynamics are indiff"erent (Di = 0). Then, [D] is positive semi-definite, and we can prove the 
global asymptotic stability of the synchronization V^^^x by using Barbalat's Lemma (see [22]). This 
special case agrees with the average consensus problem on graphs. 
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By exploiting the hierarchical structure (Theorem VIL2) of the composite variable defined in (24), we can 
show that V^^^x and A > also make qi , • • • , synchronize exponentially (i.e., V^^^ (cj^ , • • • , qj) - 
0). This straightforwardly follows the proof of Theorem III.l. 

In the case of identical spacecraft, this synchronization result also implies that the diagonal terms of 
the metric, Vf^^^[M] [1] tend to zero exponentially, thereby eliminating the couphng of the inertia term 
V^[M]V in (44). So far, we have assumed that cid{t) is identical for each spacecraft. If c^dit) were different 
for each dynamics, Sj would imply the synchronization of q^ — q^ to the difference of the desired 

trajectories, which would be useful to construct phase synchronization of spacecraft positions. Such phase 
synchronization is discussed in the subsequent section. It should be mentioned that global asymptotic 
convergence of a linear coupling control law for the nonlinear attitude dynamics is proven in Ref. 22, and 
compared with the proposed nonlinear control law in Sec. VI. 



C. Concurrent Synchronization of a Large Formation 

Consider a complex spacecraft network that consists of multiple identical or heterogeneous groups as shown 
in Figures 1(b) and 3(c). We exploit the fact that there exist two different time scales of the proposed 
synchronization tracking control law (Di and D2). In other words, we can exploit a desired reference 
trajectory c[d{t) to create multiple combinations of different dynamics groups. For instance. Figure 3(c) 
shows both the diffusive coupling structure proposed by the tracking control law (solid lines) and the reference 
directional input (dashed lines). Each spacecraft on the same network initially receives a different desired 
trajectory input from the adjacent spacecraft of the different (inner) network. Once the inner network 
synchronizes, the outer network also ends up receiving the same desired trajectory to follow, while they 
interact to synchronize exponentially fast. Accordingly, we can achieve concurrent synchronization^^ between 
multiple network groups. The proof straightforwardly follows Theorem VII. 2. Readers are referred to Ref. 22 
for further discussions. 



IV. Decentralized Nonlinear Control for Relative Translational Dynamics 

We introduce a new method for the phase synchronization of multiple spacecraft positions that follow a 
time-varying circular or spiral trajectory with equal spacing between spacecraft, as shown in Figure 4. We 
also illustrate that the rotational and translational dynamics can be combined to enable coupled rotational 
maneuvers, which can realize a spiral trajectory. The UV coverage of an interferometer defines the size 
and density of spatial points on the two-dimensional Modulation Transfer Function (MTF) of an image. ^'^ 
By spiraling out, a stellar interferometer can construct a virtually large telescope, thereby accomplishing a 
suitable UV coverage, which in turn results in a finer angular resolution and higher contrast ratio. The kinds 
of stellar objects that will be observed will also contain information at many spatial frequencies, thus the 
ability of stellar interferometers to observe at multiple baselines will be key. 



A. Rotational Phase Synchronization for Spacecraft Positions 

This section extends the phase synchronization of a single variable from Ref. 41 to Lagrangian systems 
of multiple degrees-of-freedom on a three-dimensional plane with a hierarchical control law. A new proof 
for phase synchronization is presented such that the angular transformation terms in the Laplacian matrix 
disappear. By combining the trajectory control with the synchronization coupling, we can achieve much 
more efficient and robust performance that is essential to precision formation flight of spacecraft. 

For the relative translational dynamics presented in (21), the following decentralized tracking control 
law with two-way-ring symmetry is proposed for the i-th spacecraft in the circular network comprised of p 
identical spacecraft {1 < i < p): 



Mrr,^ + Cr,,, + D(r,)r, + g(r,) - ^is^ + ^2T(^)sn + k2T^{0)s^^ 



(45) 



where ki and k2 are positive constants, and the rotation matrix T(^) G SO{2) is defined as 



T(^) 



cos i9 
1 
sin6> 



— smf 


cosO 



(46) 
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Note that T^((9) = T'^iO) = T(-6>) and T{pO) = I since = 27i/p. 




Figure 4. Geometry of a circular rotating trajectory. The x-z axes denote in Figure 2. 

For two-spacecraft networks, the last couphng term with the {i + l)-th member is not used. Also, the 
shifted reference velocity vector Yr^i and the composite variable s^°^ are defined as follows: 

rr,i = ^d,i + K^d.i - i-i), sP°' = - rr,i = - ^d,i + K^i - ^d.i) (47) 

where A is a positive constant. The time- varying desired trajectory r^^^^ is different for each spacecraft since 
its phase is shifted from its neighbors by 

= (^a(t) coscjt, yd{t), a(t)s\nuot^ =Vd 

rd,2 = {a{t) cos [uot + 6>], yd{t), a{t) sin [ujt + 6>]) = T{0)rd 

: (48) 
rd,i = (a{t) cos [cc;^ + (z - l)l9], ^^(t), a(t) sin [uot + (z - l)l9]) = T((z - l)0)rd 

Note that ^d(t) can be a time- varying function, thereby constructing a three-dimensional trajectory. In 
the case of a two-dimensional circular trajectory orthogonal to the orbital plane, the radius a{t) becomes a 
constant and yd{t) reduces to zero. 

The closed-loop dynamics for the i-th satellite, constructed from (21) and (45), becomes 

Msr + Csf + - k2T{0)s^^\ - k2T^{e)s^l\ = (49) 

By inspecting Figure 4, we can find the new flow-invariant set of s^°^ that is phase-shifted from the 
original invariant set s^°^ = 82°^ = • • • = s^^^: 

gpos ^ T^(^)sP°^ = ... = T^((p-l)^)sP°^ (50) 
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where the matrix of the orthonormal eigenvectors V is obtained from (36). 
Left-multiplying (49) with T^((i — 1)^) for each spacecraft results in 



^Ki,-K2J 



fPOS 







where 



pos] 



Mt,i 




Ct,i 




^pos 



(51) 



(52) 



Vt^((p-1)^)sP°7 



Also, Mr,i = T^((i - l)6')MT((i - 1)0) = M and Cr,i = T^((i - l)6')CT((i - 1)6') where the M and C 
matrices are defined in (21). Hence, Mr,j — 2Cr,i = — 2CT,i is skew-symmetric. Then, similar to the previous 
section, we can prove the synchronization V^„gXP°^ ^ as well as the tracking convergence s?°'^ ^ by 
constructing the virtual system and performing the spectral decomposition as follows: 



(V^[MP°'^]V) y + (V^[CP°^]V) y + [D]y = 



(53) 



Note that (53) has y = V-'^xP"^ = ([1], Vsj,„c)'^xP°^ and y = as particular solutions. Due to the skew- 
symmetry of C, (53) is contracting with [D] > 0, resulting in the tracking stability s^°^ — > as well as 

sP°« ^ T^(^)sP°" ^ T^((i - 1)^)8?°" ^ T^((p - l)^)sP°^ (54) 

From the hierarchical combination of the composite variable and the definition of r^j j, one can verify that 



h + Xvi + T((i - 1)61) (fd + Xvd) = s: 
fi+i + Ari+i + T{i9){rd + Xra) = s. 



pos 



(55) 



pos 

+ 1 



Left-multiplying the dynamic equation of in (55) with T^{{i — 1)6>), and that of r^+i with T^{{i)6) 



result in 



T^((i - l)0)ri + AT^((i - l)0)ri = u{t) 
T^(i^)r,+i + AT^(i^)r,+i = u{t) 



(56) 



where the common input verifies u{t) = — (r^^ + Ar^^) if the phase-shifted composite variables synchronize, 
that is T^{{i - l)^)sP°' = T^{iO)s^l\. 

From (56), A > ensures that T^((i — l)0)ri exponentially tends to T^(i^)r^+i since y + Ay = is 
contracting (Theorem VILl and Theorem VIL4). This concludes the proof of the phase synchronization 
of the relative translational dynamics for the case of multiple identical spacecraft on a circular or spiral 
trajectory. For heterogeneous spacecraft, we can scale the control gains ki and k2 according to the mass 
ratio, and the proof essentially remains the same. Furthermore, we can straightforwardly extend the constant 
phase shift 6 to arbitrarily different phase shifts among spacecraft (i.e., 6i ^ 6j with i ^ j). 

B. Synchronized Rotational Maneuvers 

One can easily construct the combined dynamics of both attitude rotation and relative position as follows 



■[M] 







"[C] 











[C^^^] 













Ki ,— K2 



= 



(57) 



where we used the control law (23) for the rotational attitude dynamics and the control law (45) for the 
translational dynamics. 

We can also demonstrate the synchronized rotation maneuver^^ by synchronizing the desired rotational 
rate of the attitude dynamics from (id{t) and qd(t) with the rotational rate oj of the desired circular trajectory 
Vdit) defined in (48). 
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V. Effects of Communication Delays and Model Uncertainties 



Now let us discuss the robustness properties of the proposed synchronization tracking control approach. 



A. Synchronization with Time Delays 

Based on Ref. 29, we show herein that contraction properties are conserved in time-delayed diffusion- like 
couplings for multiple Lagrangian dynamics of spacecraft. In particular, the proposed synchronization cou- 
pling control law in (23) and (45) can be proven to synchronize multiple dynamical systems as well as to 
track the common trajectory, regardless of time delays in the communication. 

Let us consider two spacecraft transmitting their attitude state information to each other via time-delayed 
transmission channels, as shown in Fig. 5. While T12 is a positive constant denoting the time delay in the 
communication from the first spacecraft to the second spacecraft, T21 denotes the delay from the second 
spacecraft to the first spacecraft. 




Figure 5. Synchronization of two identical spacecraft with transmission delays 

Similar to Ref. 29, we can modify our original Lagrangian systems consisting of two identical spacecraft 
in (28) as follows 

M(qi)si + C(qi, qi)si + (Ki - K2)si - GsiTsi = 
M(q2)s2 + C(q2, q2)s2 + (Ki - K2)s2 - Gi2ri2 = 

where G21 and G12 are constant matrices (R^^"'). 

The communication between the two dynamics occurs by transmitting intermediate "wave" variables, defined 
as^^ 

U21 = G^iSi + /C2ir2i V12 = G^iSi 

U12 = G12S2 + A:i2ri2 V21 = G12S2 

where /C21 and ki2 are strictly positive constants. Time delays of T21 and T12 result in 

Ui2(t) = Vi2(t - T12) U2i(t) = V2i(t - T21) (60) 

Notice that the original dynamics without any communication delays are contracting since Ki — K2 > 0. 
Expanding (58), using the above relationships on the wave variables, yields 

M(qi)si + C(qi, qi)si + (Ki - K2)si - -^G2i {^^2^2^ - T21) - G^iSi(t)) = 

\- (61) 
M(q2)s2 + C(q2, q2)s2 + (Ki - K2)s2 - — G12 {G^Mt - T12) - Gj^S2{t)) = 

We can verify that (61) becomes equivalent to the original two-spacecraft dynamics that can be written as 
in (28) by setting ki2 = /C21 = 1 and G12 = G21 = a/K2- Note that >/K2 is the Cholesky decomposition of 
the positive- definite symmetric matrix, K2. 

The resultant equations reflecting the time-delayed transmissions become 

M(qi)si + C(qi, qi)si + KiSi - K2S2(t - T21) = 
M(q2)s2 + C(q2, q2)s2 + K1S2 - K2Si(t - T12) = 
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which can be shown to be semi-contracting (asymptotically stable) using the following differential length 
similar to Ref 29: 

V = Sz^ (V^[M] V) Sz + 1/1,2, z = V^(sf , s^)^ (63) 

where 



Jt-Ti2 Jt-T2i 



(64) 



In conclusion, the formation flying spacecraft systems, individually contracting (exponentially converging) 
and interacting through time-delayed diffusion-like coupling are asymptotically contracting regardless of the 
values of the time delays. This shows that the proposed control law and its closed-loop system in (28) possess 
some robustness properties with respect to time delays. 

B. Effect of Bounded Disturbances 

Equation (39), in the presence of the external disturbance torque Text, can be written as 



(V-[M]V) ( ^ (vT[c]v+ [D]) ( ] ^ 

\ ^ sync J \ ^ sync J \ ^ sync j 



\Text,p/ 



(65) 



which indicates that the disturbance input for the synchronization is only the difference among each distur- 
bance force/torque V^^^(rext,i"^, • • • , Text,p"^)"^- As a result, the disturbance torque that is invariant from 
spacecraft to spacecraft, does not affect the synchronization of the relative attitude, which might be of more 
importance than the performance of trajectory following. For example, stellar interferometers need precise 
control of relative attitude and distance between spacecraft that carry telescopes. 

Now we consider the bounded vanishing disturbance of the individual tracking dynamics. Due to expo- 
nential tracking convergence of the proposed scheme, the property of robustness to bounded deterministic 
disturbances can easily be determined. For example, consider the closed-loop system in (28), which is now 
subject to a vanishing perturbation^^ such that d(t, x = 0) = 0: 

[M]x + [C]x + [L^^,_kJx = d(i,x) (66) 

The perturbation term d(t, x) vanishes at the equilibrium manifold x = 0. Let us further assume that it 
satisfies the linear growth bound such that 



||d(i,x)|| <7||x||, Vt>0. (67) 



where 7 is a positive constant. 
The squared-length analysis yields 



^ (^x^[M]^x) = 2^x^[M]^x + ^x^[M],5x 

= 2^x^(-[C](5x- [L^^__Kj(5x + (5d(i,x)) +^x^[M](5x (68) 
<-25x^([L^^,_kJ-7i) 

where we used the skew-symmetric property of [M] — 2[C]. 

Hence, the closed-loop system in (66) is contracting in the presence of the bounded disturbance if 
[L^^ -K2] ^ ^ result, the trajectory tracking gain (slower gain) also determines how robust the 

closed-loop system is with respect to a bounded disturbance. 

For a nonvanishing perturbation such that 

||d(t,x)||<7||x||+A (69) 

where A is bounded, the comparison method^^ can straightforwardly be developed to derive a bound on the 
solution. Alternatively, we can follow the analysis of contraction with respect to a property of robustness. ^^'^^ 
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We assume that Pi{t) represents a desired system trajectory and P2{t) the actual system trajectory in a 
disturbed flow field given in (66). Also, consider the distance R{t) between two trajectories Pi(t) and P2{t) 
such that ^ ^ 

R{t) = [ ' ||,5z|| = / ' ||0(x)^x|| (70) 

JPl JPi 

where 0(x)-^0(x) = [M] (see the Appendix). 

Then, by combining the results from (68), (69), and Ref. 26, we can show that any trajectory converges 
exponentially to a ball of radius of R around the desired trajectory such that 

R{t) < sup,,,||0(x)-^A||/A^«, (71) 

where the contraction rate Xmax, in the context of contraction theory, is defined with respect to Sz = 0(x)(5x, 
and is given by [M]~^ (l-'^Ki -K2] ~ from (68). It should be emphasized that the exponential stability of 
the closed-loop system facilitates such a perturbation analysis, which showcases another benefit of contraction 
analysis. In contrast, the proof of robustness with asymptotic convergence is more involved. 

VI. Extensions and Examples 

Let us examine the effectiveness of the proposed control law in a few examples. 

A. Synchronization with Partial Degrees-of-Freedom Coupling 

In this section, we consider multiple spacecraft with partially coupled variables. For instance, we can assume 
that only the first and third MRP variables {qi and gs) are coupled in a two-spacecraft system such that 



Ti = M(qi)qi^ + C(qi,qi)qi^ + g(qi) - KiSi +K2 fg^ ^3) + K2A ^3) 

\ ^ ^ (72) 

T2 = M(q2)q2r + C(q2,q2)q2r + g(q2) -K1S2 + K2 (^g^ ^3) + K2A (^g^ §3) 



Nonetheless, Theorems III.l and III. 2 are true with diagonal matrices, Ki, K2 and A, which can be 
verified by writing the closed-loop system as: 

M(qi)si + C(qi, qi)si + (Ki + K2P)si = u{t) 

M(q2)s2 + C(q2, q2)s2 + (Ki + K2P)s2 = u{t) (73) 

u(t) =K2P(S1 + S2) 

where P = diag(l, 0, 1). 

We can prove that Theorems III.l and III. 2 still hold by showing that Ki + K2P and Ki — K2P are still 
uniformly positive definite, enabling exponential synchronization and exponential convergence to the desired 
trajectory, respectively. Hence, we did not break any assumptions in the proof of Theorem III. 2. 

B. Attitude Synchronization of Two Spacecraft 

The proposed control law in (23) is simulated for two identical spacecraft, as shown in Figure 6(a). The 
spacecraft inertial matrix is 

150 -100" 
270 
-100 300 



Js/c 



[kgm^] (74) 



The control gains are defined as Ki = 3001, K2 = 1001, and A = 201. The reference trajectories are defined 
as qic^ = 0.3 sin (27r(0.01)t), q2d = 0.2 sin (27r(0.02)t + 7r/6), and c^sd = 0- The first spacecraft is initially 
at (0.05, —0.1,0)-^ rad, with zero angular rates, while all the initial conditions for the second spacecraft are 
zero. 

The synchronization gain of si and S2 corresponds to Ki + K2 = 4001, which is larger than the tracking 
convergence gain Ki — K2 = 2001. As a result, we can see in Figure 6 that the first and second spacecraft 
exponentially synchronize first. Then, they exponentially converge together to the desired trajectory. 
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Synchronization of Two Spacecraft, Proposed Approach 
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(a) Proposed control law 
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(c) Control effort: proposed control 



Synchronization of Two Spacecraft, PD Control 
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(b) Simple PD diffusive coupling 



Control Effort, PD Control 
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(d) Control effort: PD control 



Figure 6. Synchronization of the attitude dynamics of two identical spacecraft 



In order to compare the effectiveness of the exponential tracking, a simple Proportional and Derivative 
(PD) diffusive coupling, introduced in Ref. 4, 22, is simulated for the comparison purpose, as shown in 
Figure 6(b). The control law for two spacecraft can be given as 

n = -Ki(qi + Aqi) + K2(q2 + Aq2) 
r2 = -Ki(q2 + Aq2) + K2(qi + Aqi) 

whose global asymptotic stability with respect to a constant reference input is proven in Ref. 4, 22. For 
a fair comparison, we selected the PD gains in (75) as Ki = 1000, K2 = 300, and A = 0.3 such that 
the level of control efforts is comparable to that of the nonlinear control approach in (23). As shown in 
Figure 6(b), the PD coupling control law is not effective in following a time- varying trajectory. This is 
because simple linear control cannot be expected to handle the dynamic demands of efficiently following 
time-varying trajectories. In contrast, the proposed nonlinear approach shown in Fig. 6(a) exponentially 
synchronizes the attitude states, while the tracking errors tend exponentially to zero. Specifically, achieving 
exponential convergence ensures more effective tracking performance than asymptotic convergence by linear 
PD control. The proposed nonlinear control law with the given initial condition introduces a large control 
effort initially. However, the results do not change, even if the torque actuators saturate at ± 6 Nm. 
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C. Four Spacecraft Example 

In order to illustrate the effectiveness of the proposed approach for a larger formation, a result of simulation 
for four non-identical spacecraft is presented in Fig. 7. The spacecraft inertia matrices are 



150 
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100 
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in [kgm^], 


(76) 
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respectively. The desired trajectories are the same as in the previous section, and the control gains remain 
the same as well, such that the tracking convergence gain is Ki — 2K2 = 1001 > 0. Figure 7 shows that the 
four spacecraft synchronize themselves while following the desired trajectory together. 



D. Position and Attitude Synchronization of Two Spacecraft Under J2 Effect 

A full 6 DOF simulation of a two identical satellite formation using the proposed phase synchronization 
control law (45) and with the attitude tracking (2) was prepared. Full nonlinear attitude and translational 
dynamics including J2 effect were simulated for each satellite in the formation. The mass of each satellite 
is assumed to be 500 kg. The spacecraft inertia matrix is given by (74). The attitude control gains were 
defined as Ki = 301, K2 = 201 and A = 201. The attitude reference trajectories were defined to be 
Qid = 0.3sin(27r x O.OOlt), q2d = 0.2sin(27r x 0.002t + 7r/6), and qsd = 0, which are similar (but at a lower 
frequency) to the previous example. 

For the translational control the gains were defined to be same for both the satellites as /ci = 10, /c2 = 5 
and A = 2. The desired trajectory for the formation was defined so that the satellites would follow a spiral 
trajectory such that both satellites maintain a phase angle of 180 deg. The desired spiral trajectory was 
defined as r^ = [a{t) sinc^t, 0, a{t) cos cot]^ from (48), where a{t) = 5 + O.OOOlt and co = 27r/500 (see the actual 
plot of this spiral trajectory generated from the actual simulation in Fig. 8). The formation was placed in a 
polar circular orbit at an altitude of 500 km and it was subject to J2 perturbations in the nonlinear model. 
Figure 9 shows the robust synchronization performance in the presence of non-identical disturbance such as 
the differential J2 effect. It should be noted that this scheme requires information from only the adjacent 
spacecraft, thus this scheme is suitable for implementing a distributed control architecture. 

VII. Conclusions 

We have introduced the new unified synchronization framework that integrates both the exponential 
tracking of a demanding time- varying trajectory and the synchronization of spacecraft motions either through 
local coupling feedback or a relative sensing metrology system. The new decentralized control laws, developed 
by utilizing the Lagrangian formulations of the translational dynamics and attitude dynamics of spacecraft, 
enable coupled rotational maneuvers and phase synchronization of the attitude and position, thereby facil- 
itating a further analysis on stability, adaptation, and robustness. In order to rigorously address the three 
main areas of research for the realization of future spacecraft formation flying missions, we have focused on 
the three research areas: (1) the exact nonlinear stability conditions of multiple spacecraft dynamics, (2) the 
reduced information networks through partial state couplings and local interactions, and (3) the properties of 
robustness with respect to uncertain models and time delays of the communication. In particular, in contrast 
with prior work which used simple single or double integrator models, the proposed method is applicable to 
highly nonlinear systems with nonlinearly coupled inertia matrices such as the attitude dynamics of space- 
craft. It should be noted that others can still come up with nonlinear control laws, similar to our approach, 
without a rigorous stability proof. Hence, one of the main contributions of this paper lies with the use of a 
new differential nonlinear stability framework called contraction theory, yielding the exact proof of nonlinear 
stability under a variety of conditions. Contraction analysis, overcoming a local result of Lyapunov's indirect 
method, yields global results based on differential stability analysis. The benefit of constructing multiple 
timescales of the closed- loop system is that exponential synchronization, with a convergence rate faster 
than that of the trajectory tracking, enables reduction of multiple dynamics into a simpler synchronized 
form, thereby simplifying the additional stability analysis. The proposed bi-directional coupling has also 
been generalized to permit partial-state coupling, thereby further reducing communication requirements. 
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Synchronization of four heterogeneous spacecraft 
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Figure 7. Synchronization of four heterogeneous spacecraft. 
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Figure 8. Three-dimensional simulation of the phase synchronization on a spiral trajectory 



10 

f 
^-10, 

10 

^ 
-10, 

10 

I 

N 

-10 



(a) Results of the simulation for two spacecraft phase syn- 
chronization 
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Figure 9. Simulation of the proposed control law for the combined translational and rotational dynamics. 
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We illustrated the effectiveness of the proposed approach by comparing with the simple Proportional and 
Derivative (PD) coupling control law, and simulating the combined rotational and translational maneuvers 
in the presence of variation in the second-degree zonal harmonic J2 of the Earth's gravitational potential. 
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APPENDIX: Contraction Theory 

We exploit partial contraction theory^^ to prove the stability of coupled nonlinear dynamics. Lyapunov's 
linearization method indicates that the local stability of the nonlinear system can be analyzed using its 
differential approximation. What is new in contraction theory is that a differential stability analysis can be 
made exact, thereby yielding global results on the nonlinear system. A brief review of the results from^^^^ is 
presented in this section. Readers are referred to these references for detailed descriptions and proofs on the 
following theorems. Note that contraction theory is a generalization of the classical Krasovskii's theorem. 

Consider a smooth nonlinear system 

x(t)=f(x(t),u(x,t),i) (77) 

where x(t) G M^, and f : x x IR+ ^ M^. A virtual displacement, Sx is defined as an infinitesimal 
displacement at a fixed time- a common supposition in the calculus of variations. 

Theorem VII. 1 For the system in (77), if there exists a uniformly positive definite metric, 

M(x,t) = 0(x,t)^0(x,t) (78) 

where is some smooth coordinate transformation of the virtual displacement, Sz = SS:s., such that the 
associated generalized Jacobian, F is uniformly negative definite, i.e., 3A > such that 



F = (&{x,t) + 0{x,t)^) 0{x,t)-' < -AI, 



(79) 



then all system trajectories converge globally to a single trajectory exponentially fast regardless of the initial 
conditions, with a global exponential convergence rate of the largest eigenvalues of the symmetric part ofF. 

Such a system is said to be contracting. The proof is given in Ref.^^ Equivalently, the system is contracting 
if 3A > such that 



M+ h- M + M^ < -2AM (80) 
\dxj ox 

It can also be shown that for a contracting autonomous system of the form x = f (x, u(x)), all trajectories 
converge to an equilibrium point exponentially fast. In essence, contraction analysis implies that stability 
of nonlinear systems can be analyzed more simply by checking the negative definiteness of a proper matrix, 
rather than finding some implicit motion integral as in Lyapunov theory. 

The following theorems are used to derive stability and synchronization of the coupled dynamics systems. 

Theorem VII. 2 Hierarchical combination^^' 

Consider two contracting systems, of possibly different dimensions and metrics, and connect them in series, 
leading to a smooth virtual dynamics of the form 




Then the combined system is contracting if F21 is bounded. 
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Theorem VII. 3 Partial contraction^^ 

Consider a nonlinear system of the form x = f(x, x, t) and assume that the auxiliary system y = f(y,x, t) 
is contracting with respect to y. If a particular solution of the auxiliary y -system verifies a specific smooth 
property, then all trajectories of the original x-system verify this property exponentially. The original system 
is said to be partially contracting. 

Theorem VII. 4 Synchronization^^ 

Consider two coupled systems. If the dynamics equations verify 

xi - f(xi,t) = X2 - f(x2,t) 

where the function f(x, t) is contracting in an input-independent metric, thenyii and ^2 will converge to each 
other exponentially, regardless of the initial conditions. Mathematically, stable concurrent synchronization 
corresponds to convergence to a flow-invariant linear subspace of the global state space.^^ 
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